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ABSTRACT 

Low-mass stars are known to experience extra-mixing in their radiative zones 
on the red-giant branch (RGB) above the bump luminosity. To determine if the 
salt-fingering transport of chemical composition driven by ^He burning is efficient 
enough to produce RGB extra-mixing, 2D numerical simulations of thermohaline 
convection for physical conditions corresponding to the RGB case have been 
carried out. We have found that the effective ratio of a salt-finger's length to 
its diameter Ocfr < 0.5 is more than ten times smaller than the value needed to 
reproduce observations (oobs ^ 7). On the other hand, using the thermohaline 
diffusion coefficient from linear stability analysis together with a = Oobs is able 
to describe the RGB extra-mixing at all metallicities so well that it is tempting 
to believe that it may represent the true mechanism. In view of these results, 
follow-up 3D numerical simulations of thermohaline convection for the RGB case 
are clearly needed. 

Subject headings: stars: abundances — stars: evolution — stars: interiors 

1. Introduction 

A main sequence (MS) star with a mass M < 1.5 Mq, e.g. the Sun, gets its energy 
from thermonuclear fusion reactions of the pp-chains. When its core has exhausted the 
hydrogen fuel, the star leaves the MS and becomes a red giant branch (RGB) star. The red 
giant has a compact electron-degenerate helium core of the size of the Earth surrounded by 
a thin hydrogen burning shell (hereafter, HBS) that provides the star with the energy via 
the transformation of H into He via the CNO-cycle. The bulk of the red giant's volume is 
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occupied by a convective envelope that extends from the surface to nearly the HBS, however 
the base of the convective envelope always stays separated from the HBS by a convectively 
stable radiative zone with a thickness of the order of one solar radius. 

It is possible to select a large sample of low-mass post-MS stars that have very similar 
masses and metallicities. The simplest way to do this is to choose normal post-MS stars be- 
longing to a same several billion year old open or globular cluster. For field stars, one needs 
well-measured distances to estimate the masses accurately enough, so that only stars with 
well-determined Hipparcos parallaxes suit this purpose. A distribution of the characteristics 
of stars from such a sample along the RGB in the Hertzsprung-Russell (HR) diagram can 
then be interpreted as an evolutionary sequence (the Vogt-Russell theorem). In particular, 
a correlation of a surface chemical composition anomaly (a deviation of atmospheric abun- 
dances from those on the MS) with the luminosity or effective temperature can be considered 
as an evolutionary pattern produced by nuclear reactions and ongoing mixing. 

Spectroscopic observations and their theoretical interpre tations have firmly established 



two distinct mixing episodes in low-mass RGB stars (e.g., ICharbonnell Il994l ). First, on 
the lower RGB, the convective envelope quickly grows in mass and, as a result, its base 
penetrates the layers in which the partial burning of H via the pp-chains resulted in an 
accumulation of large amounts of ^He, and where Li was mostly destroyed. It also reaches 
the layers in which the non-equilibrium operation of the CN-branch of the CNO-cycle took 
place. This should lead to sharp changes of the star's atmospheric chemical composition, 
such as a strong decrease of its surface Li abundance, a noticeable reduction of the ^^C/^^C 
isotopic ratio, a modest decrease and increase of the C and N abundances respectively, and 
a considerable enrichment of the convective envelope in ^He. These changes continue until 
the convective envelope attains its maximum mass after which its base starts to recede in 
front of the HBS advancing in mass. This mixing episode is called the first dredge-up. All 
of the above predicted changes of the star's surface chemical composition, except those of 
^He, have often been observed on the lower RGB. 

Second, on the upper RGB, above the luminosity at which differential luminosity func- 
tions of globular-cluster RGB stars have clearly visible bumps, observations reveal that the 
same abundance changes that were occurring during the first dredge-up resume, and that 
the changes are even larger this time. The bump luminosity corresponds to the evolutionary 
phase during which the HBS passes through and erases the chemical composition disconti- 
nuity imprinted by the convective envelope at the end of the first dredge-up. When the HBS 
crosses the discontinuity it finds itself in a region with a slightly increased H abundance; con- 
sequently, the red giant readjusts its structure appropriately, causing it to move to slightly 
lower luminosities before continuing up the RGB. As a result of spending more time in the 
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narrow range of luminosity where this readjustment occurs, there is a pile-up of stars around 
the bump luminosity. Given that the base of convective envelope recedes through a region 
of homogeneous chemical composition, the observed evolutionary variations of the surface 
abundances on the upper RGB can be understood only if there is some extra-mixing process 
of non-convective origin that connects the convective envelope with the HBS. The physical 
nature of this extra-mixing still remains elusive. For a long time, it ha s been thought that 



the RGB extra-mixing i s related to rotation- i nduce d instabilities (e.g., ISweigart fc Mengel 



I979I: ICharbonnell Il995l iDenissenkov fc ToutI I2OOOI ). However, the work of iPalacios et al 

casts doubt on these models 



It was not u ntil recently that the most promising me chani sm of the RGB extra - mixiri g 
was proposed by lEggleton. Dearborn, fc Lattanzid (120061 ) and ICharbonnel fc ZahnI ( 120071 ). 
When the HBS is approached from the base of convective envelope, the first active ther- 
monuclear reaction involving abundant nuclei that is encountered is ^He(^He,2p)^He. It has 
a unique property of lowering the mean molecular weight locally, by A/x ~ /i^AX(^He)/6, 
where AX(^He) is a mass fraction of ^He consumed in the reaction. As a consequence, the 
density is reduced by the same amount (assuming the ideal gas law), which makes the local 
material lighter than that in its immediate surroundings, therefore the ^He-depleted mate- 
rial will tei id to rise. M ixing driven by ^He burning in s tars was discussed already in the 



1970s (e.g.. lUlrichlll972l ). however it was the publication by lEggleton. Dearborn, fc Lattanzio 



(120061 ) that brought this mechanism to the attention of other researchers in the context of 
the RGB extra-mixing. 



Charbonnel fc ZahnI (120071 ) have given a correct physical explanation of the mixing pro- 
cess to which the ^He burning in the HBS should lead. This is the salt-fingering or thermo- 
haline convection that has been studied in detail, both experimentally and theoreticall y, in 
oceanograph y for many years (e.g., see reviews by lRuddick fc Kerr II2003I . ISchmittll2003l . and 



Kunze 



2OO3I ). It is caused by a double- diffusive instability that occurs in a situation when 



a stabilizing agent (heat) diffuses away faster than a destabilizing agent (salt and /i in the 
oceanic and RGB cases, respectively). In the ocean, the double-diffusive instability usually 
develops close to the surface where warm salty water finds itself overlying cold fresh water. 
A blob of the surface water will tend to sink because its higher salinity (a concentration of 
salt) makes it denser than the surrounding deeper water while its temperature remains close 
to the ambient one thanks to the faster diffusion of heat. Similarly, a blob at depth will tend 
to rise. In the oceanic case, the usual outcome of this instability are vertically elongated salt 
fingers containing sinking and rising water of different salinity. 

Given that the local reduction of n by the ^He burning is very small, AfJ^/fJ^ ~ —10^^, 
it can become visible only in the background of homogeneous chemical composition. This 
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happens precisely at the bump luminosity. Another feature that makes the ^He-driven 
thermohaline convection the most promising mechanism for the RGB extra-mixing is its 
dependence on just one parameter, the envelope ^He abundance after the first dredge-up. 
This means that stars with similar masses and metallicities should demonstrate comparable 
chemical compo sition anomalies on the upper RGB, the patt ern that seems to be observed in 
real stars (e.g., iGratton et al.ll2000l : ISmith fc Martelll 120031 ). Before the RGB extra- mixing 
was introduced into the standard stellar evolution theory, the latter had been in conflict with 
the seemingly observed constancy of the interstellar ^He abundance since the Big Bang. The 
problem is that low-mass stars produce a lot of ^He on the MS that is dredged up to the 
surface on the RGB and then deposited in the interstellar medium as a result of mass-loss. 
The most plausible solution of this problem is the RGB extra-mixing that strongly decreases 
the envelope abundance of ^He by circulating it through the ^He burning layers of the HBS. 
Therefore, if the RGB extra-mixing is indeed driven by the ^He burning (thermohaline 
convection) then the cosmological ^He problem has a beautiful solution: it is ^He itself 
abundantly produced in low-mass MS stars that takes care of its own destruction in the 
same stars on the RGB, so that the net balance of '^He from low-mass stars in the interstellar 
medium is close to zero. 



Charbonnel fc ZahnI (120071 ) have mo deled thermoh aline mixing using a diffusion coef- 
ficient obtained from a linear theory by lUlrichI (jl972l ). Unfortunately, the linear theory 
does not give a reliable estimate of the maximum length of salt fingers relative to their 
diameter, i.e., the finger aspect ratio a = l/d, the square of which enters the diffusion coeffi- 
cient. This leads to a large uncertainty in the theory leaving it basically semi-empirical, like 
ot her up-to-date theories and models of the RGB extra-mixing. From the results presented 
by I Charbonnel fc Zahru (120071 ) , it follows that the surface abundance patterns in low-mass 
RGB stars can be reproduced the oretically only if a 3> 1. To support the large finger aspect 



ratio, I Charbonnel fc ZahnI (120071 ) referred to the oceanic case where long salt fingers were 
observed experimentally. However, it is not legitimate to directly compare the oceanic and 
RGB cases because they correspond to very different fluid flows sets. In particular, the ratio 
of viscosity to thermal diffusivity (the Prandtl number Pr) is close to ten for the oceanic 
case, whereas Pr < 10~^ for the RGB case. It turns out that only direct numerical sim- 
ulations with non-linear interactions and other relevant effects taken into account can give 
better insight into the ^He-driven thermohaline mixing in the low-mass RGB stars. It is this 
task that is addressed in this research for the first time. 

The paper is organized as follows. Section H] summarizes the main results from the linear 
theory that are relevant for our further discussion. Section IHlpresents and analyzes the results 
of our 2D numerical simulations of thermohaline convection. It is followed by Section |U in 
which we employ a diffusion coefficient from the linear theory to model, as precisely as 
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possible, the changes of chemical composition incurred by the RGB thermohaline mixing. 
Discussion and conclusions are provided in the final Section |5l Wherever it is appropriate, 
we make a comparison of the oceanic and RGB cases. 



2. Relevant Results from the Linear Salt-Fingering Theory 

In this as well as in the next section, our analysis begins with the Boussinesq equations 
that describe motion in a nearly incompressible stratified viscous fluid 

|^ + (v,V)v = ^g + z.VV, (1) 
dt po 

dT 

— + (v,V)T = A^tV^T, (2) 
|i + (v,V)^ = ksV'S, (3) 

where v is the velocity, po is a constant reference density, 5p is a deviation of the local density 
from Po, g is the gravitational acceleration, v is the viscosity, T is the temperature, 5" is the 
salinity, fc^ and ks are the thermal and haline diffusivities. Note that, in the RGB case, S 
should be replaced by p, while equation ([3]) is still valid as long as | Ap|/p ^ 1. We have also 
neglected source terms in equations ((21) and ([3]) that could be related to nuclear reactions in 
the RGB case. 

In the Boussinesq approximation, which is a reasonable one even for a compressible 
fluid provided that its motion is studied on length scales much less than the density height 
scale and velocities remain much less than the speed of sound, the density variation is 
taken into account only in the buoyancy force term fb = (5p/po)g. Assuming that the 
relative deviations of p, T, and 5* from their reference values in the initial unperturbed state 
(po,To, 5*0) are small, we use a linearized equation of state 

5p = -a6T + /36S, (4) 

where 

1 f din p\ 1/(9 In p 

a = — — — ; — , and p 



T\d\nTjp' S\d\nS, 

are the coefficients of thermal expansion and haline contraction. For the ideal gas law, which 
provides a good approximation to the equation of state in the radiative zone of a low-mass 
RGB star, we simply have a = 1/T and /3 = 1/p. 

We use the Cartesian coordinate system (x, y, z) oriented so that its vertical axis z 
has a direction opposite to that of the gravitational acceleration and the x-axis is located 
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in a star's meridional plane. Let w denote the velocity's vertical component. Linearizing 
equations ([I]-|3]), we arrive at the following system of linear PDEs for the vertical velocity 
component and variations of temperature and salinity: 

— = g{aST-l35S) + uV^w, (5) 

^ - + *.V=^S. ,7) 

Given that we seek a solution representing vertically elongated structures with a large ratio 
of the vertical to horizontal length scales, we can neglect horizontal velocity components and 
we can also assume that ~ d'^/d'^x + d'^/d'^y. 



Following iKunzd (119871 ). we will take into consideration the influence of salt-fingering 
on the ambient temperature and salinity gradients that appear on the right-hand sides of 
equations (E]) and (^^. This influence is determined by the ratios 6t = —26T/AT and 6s = 
—26S/AS, in which AT = (dTo/dz)! and AS = {dSo/dz)l are differences in temperature 
and salinity between two points separated by the finger length / in the vertical direction in the 
initial unperturbed state. The factor of two in the above relations takes account of the fact 
that every sinking finger with excesses of S and T has a neighbouring rising finger possessing 
deficiencies of S and T of same magnitudes. The perturbed gradients are expressed via the 
new independent variables 

dT dTn ,^ . , 



dz dz 

where < 6s) < 1- With regard to the temperature gradient, it is important to note 
a difference between the oceanic and RGB cases. In the RGB case, one should subtract 
the adiabatic temperature gradient {dT jdz)^^ from both dT /dz and dT^/dz in equation ^ 
because a temperature gradient is always negative in stars and it remains stable only as long 
as its absolute magnitude is less than that of the adiabatic gradient. The maximum effect 
that a salt-fingering heat flux can produce on dT /dz is to render it adiabatic {5t = 1)- In 
the oceanic case, any positive temperature gradient stabilizes the density stratification. 

Finally, expressing 5T and 5S by 5t and 5s and substituting them into the system ([5]-[7]) 
gives 



^-^V^^ = \9P^-I^{5s-R,5tI (10) 
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^-krV'Sr = jwil-6T), (11) 

^-ksV'6s = jw{l-Ss), (12) 

where Rp = aAT/ PAS is a parameter known as the density ratio in oceanography. In 
stellar physics, it corresponds to Rp = (V — Vad)/V^, where V = (dlnTo/dlnP), Vad = 
(9 In T/c? In P)ad, and V^^ = {d\nfio/d\nP). Solutions of the last equations are sought in the 
usual form of w, 6t, Ss oc exp(crt) x exp[i{kxX + kyy)] taking into account that the vertical 
velocity can be approximated as w ^ crl/2, where a is the growth rate of salt fingers, whereas 
kx and ky are their horizontal wave numbers. After some simple algebra, the three equations 
are reduced to the following third-order dispersion relationship: 



Aa^ + 2{kT + ks + 2u)k^a^ + {[krks + 2v{kT + ks)]k^ + 

2gP (Rp -l)}a-gP-^ (k^ - Rpks)^ = 0, (13) 



where k = ^k"^ + ky is related to the salt-finger diameter d = 2t[ /k. 



Positive real roots of equation ( 1131) . corresponding to exponentially growing solutions 
(the double-diffusive instability), are plotted in Fig. [T]for the oceanic (upper panel) and RGB 
(lower panel) cases. Note that, from here on, we will be using more customary notations for 
the stellar structure parameters related to salt-fingering, whereas those introduced so far will 
be reserved for the oceanic case. In particular, the thermal and haline diffusivities become 
the radiative and molecular diffusivities in the RGB case 

kr ^ K 



?>xCpp^ ' 

7^5/2 

ks z^moi = 1.84 X 10"^^(1 + 7X) (cm^ s'M, 

P 

where a is the radiation constant, c is the speed of light, x is the Rosseland mean opacity, Cp 
is the specific heat at constant pressure, and X is the hydrogen mass fraction. The viscosity 
in the vicinity of the HBS consists of comparable parts of z/moi and the radiative viscosity 

4ar4 

= T^mol + t'rad, where Z/rad 



2 ■ 



15cxp' 

For the reader's convenience. Table [T] shows the correspondence between similar parameters 
in the oceanic and RGB cases as well as their characteristic values. For the oceanic case, 



we h ave used the data compiled by iKunzd (119871 ) (also, see Table 1 in the review by iKunze 



20031 ) . For the RGB case, the data are taken from our bump luminosity model that has a 



mass M = 0.83 Mq and heavy-element mass fraction Z = 0.0005. They refer to a point 
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in the ^He-burning region of the HBS shell at which is negative and has its maximum 
absolute value (for model details, see Section 3.3). 

From the lower panel of Fig. [H it is seen that Kk"^ ^ 2a and Kk'^ ^ 2vk'^ ~ Vmo\k'^ in 
the RGB case. This allows one to approximate the cubic equation (!T3|) by the quadratic 



K{2u + u^,{)k^ + 2gl3^{Rp 



a 



0. 



(14) 



For short horizontal wave numbers, such that 

^4 ^^|fo (^^^P^^ andA:^«2,/^° ^"^^ 

a positive solution of ( 11^ is simplified to 



dz K(2i> + Z/mol) 



cr 



I {K - RpU^oi) 



11 



k' 



2r 



d?' 



The last expression describes very well the decline of the finger's growth rate with an increase 
of its diameter ai c/max = d{a-a,p,^). It can be used to estimate a thermohaline diffusion 
wl 



coefficient 



'maxy 

~ \al'^- This leads to the well-known results 



1 {K - RpVraox) , 2,2 



4 



Rn 



kH' = C-. 



V,. 



rad 



Vad - V, 



^rao\ Vrad ^ad 



(15) 



where C = Cxunze = vr^ and 



rad 



(9 In To 
d\nP 



IQ-nGac Mr 



is the radiative temperature gradient, and Lr being the mass and luminosity at the radius 
r in the star (assuming spherical symmetry). The first of the expressions ( |T5l) indicates that 
the double- diffusive instability develops only when 1 < Rp < 1/r, where r = h'rna] /K = 
ks/kr is the inverse Lewis number. This result is known since the pioneering work of IStern 



( 19601). Its discussion in regard to the RGB thermohaline m ixing has bee n presented by 
Denissenkov fc Pinsonneaultl (l2008bl ) (also, see IVauclairl 120041 ) . lUlrichI ( 119721 ). who neglected 
the influence of salt-fingering on the ambient T- and /i-gradients, had obtained an expression 
for Dp with the constant C = Cuirich = 87r^/3. He did not retain the ir npeding factor in the 
parent h eses either. It is U l rich's formula that has been employed by ICharbonnel fc Zahn 
( 120071 ). iKippenhahn et al.l ( ll980l ) used a different physical approach to arrive at a similar 
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expression for with C = Cxipp = 12. In addition, they came to the conclusion that, 
because of their strong interactions with the surrounding medium, the perturbed blobs of 
fluid did not have a chance to become salt-fingers, so that their aspect ratio should be as 
small as a = l/d ^ 0.5. 

It is useful to estimate the diameter of the fastest growing fingers (imax as a function 
of other parameters. This can be done by analyzing the dependence of the positive root of 
equation fll4p on the wave number k. It turns out that a sufficiently good approximation for 
both the oceanic and RGB cases is 



where is the square of the unperturbed T-component of the Brunt- Vaisala (buoyancy) 
frequency. For the RGB case, we have N"^ = (Vad — ^ra.d)g/ Hp, Hp being the pressure scale 
height. We find that ^ 1.5 x 10~^s~^ at the same location in our RGB model for which 
the other data are listed in Table [H Its substitution into (fT6!) together with the values of z/ 
and kp = K gives log^g c?max ~ 4.4 which is very close to the position of the maximum on 
the solid curve in the lower panel of Fig. [H 

Solid curves in Fig. [1] demonstrate that the diameter and velocity of the fastest grow- 
ing salt fingers in the oceanic {Rp = 1.6) and RGB cases are cimax'^" ~ 3 cm, wSax'^" ~ 



0.003 a^^x""" cm s"\ and ^ 0.3 km, w^^^ ^ 2a^^^cms"\ where a^^^'' and are 



the (unknown) finger aspect ratios for the corresponding cases. The major uncertainty in the 
linear theory is the parameter a. Its value can be estimated only by numerical simulations 
that directly solve the original non-linear equations ([T]-[3]). The double- diffusive instability 
is a primary instability. It is responsible for the initial growth of salt fingers. What will 
happen later on and, in particular, how far the fiuid blob will travel vertically, forming a salt 
finger before its trajectory is bent or the finger gets destroyed by its interactions with other 
fingers and surrounding medium, is determined by secondary instabilities. It is this problem 
that we are going to address in the next section. 




(16) 



3. 



2D Numerical Simulations of Thermohaline Convection 



3.1. Basic Equations and a Method of Their Solution 



In the Boussinesq approximation, the continuity equation is simplified to (V, v) = 0. It 
can be automatically satisfied by representing the velocity vector with a stream function ip, 
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such that 

v^(.»)^(-||,f), (IT) 

where u is the velocity (horizontal) x-component. After the substitution of ( !T7|) into the 
original system ([T]-[3]), we obtain a new system of Boussinesq equations that contain only 
scalar functions and their derivatives 

= -J(^,V',P) + Pr^(T' - ^) + PrVSp, (18) 

f - -i(^,r')-f .v^r', (19) 

f - -mS')-^^rV^S', (20) 

where J {a, b) = — is the Jacobian. Unlike the quantities 5T(x, y, t) and 6S{x, y, t) 
in equations (l5]-[7l), the variances T' and S' describe deviations of temperature and sahnity 
from To{z) and So{z) not only in the horizontal plane but also along the radius, so that 

T{x, z, t) = To(0) + (^^^ z + T'ix, z, t), 

dSc 
dz 



Six, z, t) = 5o(0) +[^] z + S\x, z, t), 



where To(0) and 5*0(0) are constants measuring the temperature and salinity in the initial 
unperturbed state at an arbitrarily chosen vertical position that corresponds to z = 0. 
Equations (fTS] - [2(I]) have been non-dimensionalized using d and /kx as the length and time 
units, where d = (iinax/2vr, and dividing T and S by (dTo/dz) d and (dSo/dz) d, respectively. 



We solve the system of equations f[T8]- [20|) numerically for a 2D doubly-periodic Carte- 
sian domain by employing a computer code kindly provided by Bill Merryfield. According 
to hi m, the code uses a Fourier collocation method with dealiasing that follows the 2/3 
rule ( ICanuto et al.lllQSq ). Integration in time is via the leapfrog method, with the time- 
splitting instability damped by app lying a Robert filter with the parameter 0.002 (e.g., 
Whitfield. Holloway. fc Holyerlll989l ). Dissipation terms are represented by exponential in- 
tegration factors. The linear portion of the code was tested by co mparing the e volution of 
small initial perturbations to predictions of linear stability theory ( 1Schmittlll979l ). To check 
the non-linear portion, the code was tested for conservation of temperature variance, salin- 
ity variance, energy and enstrophy after the dissipation, buoyancy, and background gradient 
terms had been removed. Initial S' is specified by selecting Fourier coefficients from a bi- 
Gaussian distribution and scaling by where k is the magnitude of wave vector k. The 
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variance of initial non-dimensional S" is normalized to 0.1. Initial values of T' and ip are set 
to zero. 



3.2. The Oceanic Case 

For test purposes, we have first reproduced one of the results presented by Merryfiled & 
Grinder (2000) in their unpublished paper (personal communication from Bill Merryfield). 
It corresponds to the oceanic case with the density ratio Rp = 1.6 (Table [1]). As output, the 
code directly gives the effective fingering salt and heat diffusivities, Ds (denoted as in 
the RGB case) and Dt, normalized by {K). This is achieved by a procedure that first 
determines the fingering salt (/z) and heat fluxes by spatially and temporally averaging the 
products wS and wT and then dividing them by dSo/dz (d^o/dz) and OTq/Oz, respectively, 
after the fingers have attained statistical equilibrium. The horizontal size of the solution 
domain is initially chosen to be wide enough, = Lz = 16 (imax, so that the averaging does 
not produce large-amplitude fluctuations. The blue and red solid curves in the lower panel 
of Fig. [2] show transitions of the ratios Ds/kx and Dx/kj- to their equilibrium values. The 
upper panel gives a snapshot of salt fingers in the statistical equilibrium. This solution has 
been obtained with a 1024 x 1024 resolution for the salinity, while the code always uses a 
2x2 times lower resolution for T and ip. The dashed blue line in the lower panel plots the 
constant Ds/kx = 27r^a^/i?p, where a = 3. It approximates the linear thermohaline diffusion 
coefficient (ITS!) for C = 2CKunze = 27r^ = 19.7, given that r ^ 1 (Table [1]). Here, as well as in 
the next section, we will compare our numerical results for Ds and with those predicted 
by the linear theory (eq. [15]) using C = 2CKunze- With this choice, C lies approximately 
in the middle between Cxipp = 12 and Cuirich = Svr^/S = 26.3. The comparison gives an 
effective value of the finger aspect ratio that, after being substituted in (ITB]) . leads to a value 
of the linear diffusion coefficient of the same order of magnitude as the one derived from 
our numerical simulations. In particular, we have obtained a ^ 3 for the test case. This 
estimate, as well as the salinity patterns in the upper panel of Fig. |2l confirm that we really 
deal with (i.e., reproduce numerically) vertically elongated structures (fingers) in the oceanic 
case. 

Note that, although not perfect, the 2D numerical simulations of salt-fingering in the 
ocean have succeeded in producing both the salt and heat fluxes compatible with those 



measured by lSt. Laurent fc Schmit 



(e.g.. Fig. 9 in the review by 



Kunze 



20031 ) 



(119991) in the North Atlantic tracer release experiment 
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3.3. The RGB Case 



We have employed two stellar evolution codes to compute bump luminosity models 
of low-mass RGB stars: our own code, the most recent version of which is described by 



Denissenkov et al.l ( 120101 ). and the MESA code available at http : //mesa, sourceforge .net , 
The combinations of mass, heavy-element and helium mass fractions for which the models 
have been computed with our code are (M/Mq, Z, Y) = (0.83, 0.0005, 0.24), (0.83, 0.001, 0.24), 
and (0.83, 0.002, 0.24), while the models produced with the MESA code have {M/Mq, Z, Y) = 
(0.85,0.0001,0.24), (0.83,0.001,0.24), and (1.5,0.0188,0.27). The MESA lowest-metallicity 
model reproduces conditions at which the RGB e xtra-mixing; is thought to be mo st efficient, 
according to both observations and theory (e.g., iMartell. Smith, fc Brileyl l2008l ). The pa- 



rameters of the common model, M = 0.83 and Z = 0.001, are close to those of the 
field stars with known Hipparcos parallaxes for which e volutionary abundance changes have 
been found on the upper RGB by lGratton et al.l ( 20001 ). Finally, the solar-metallicity model 
has a mass and comp osition typical for the so-called Li-r ich giants that are located close to 
the bump luminosity ( Charbonnel &: Balachandran 20001) and in which Li is believe d to be 
produced in a large amount by enhanced extra- mixing (jPenissenkov fc Herwid 120041 ). 



Sharp declines of both the Li abundance and ca rbon isotopic ratio at the burnp lu- 



minosity in the majority of low-mass RGB stars (e.g., iGratton et al.ll2000l : IShetrond 12003 



D'Orazi fc Marinol 120101 ) demonstrate that the RGB extra- mixing starts to operate and be- 
comes very efficient already at this evolutionary phase. Hence, if it is actually caused by the 
^He burning and its associated thermohaline convection, the latter should already be present 
in the bump luminosity models in which the HBS has erased the chemical composition dis- 
continuity left behind by the base of convective envelope at the end of the first dredge-up 
and, as a result, the fi depression produced by the ^He burning is now a prominent feature 
on the otherwise uniform /i profile outside the HBS. 

MESA is a state-of-the-art stellar evolution code that even allows non-stop computations 
through the helium core fiash toward the end of the asymptotic giant branch (AGB) evolution 
for low-mass stars. We will use the MESA models as a background in our post-processing 
ID simulations of the RGB extra-mixing. A comparison shows that the common models 
in the two sets of computations have very similar structures. Curves in the lower panel 
of Fig. E] show the density ratio profiles in the vicinity of the HBS in our models. Their 
minima are close to the value of Rp = 1700 that is used in our 2D numerical simulations 
of the RGB thermohaline convection. They correspond to the maximum negative values 
of Vfj, ~ — 10~^ that are reached a short distance outside of the fi depression fioor, where 
= and Rp = oo. The approximate values of other relevant parameters at the location 
of the minimum Rp are listed in Table [TJ They do not vary appreciably between the models 
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with different values of Z. 



Denissenkov fc Pinsonneaultl ( l2008bl ) have shown that evolutionary changes (those which 



vary with luminosity) of the surface abundances of Li, C, N, and of the i2q^i 3q y&Xao above 



the bu mp luminosity in the field RGB stars with Z ^ 0.002 from the study of iGratton et al. 



(120001 ) can be reproduced theoretically if extra-mixing in their radiative zones is modeled 
using a diffusion coefficient D^^^ ^ 0.02K. This gives an order-of-magnitude estimate of 
the rate of extra-mixing that has to be provided by its correct physical mechanism. The 
value of D^/K = D^^^/ K = 0.02 is plotted in Fig. [3] (green line). It is to be compared 
with the red curve in the same figure that has been computed using the same method and 
resolution which were used to prepare the blue curve in Fig. [2] for the oceanic case but now 
for the set of parameters corresponding to the RGB case (Table [T]). To find out how the 
numerical value of depends on the density ratio, we have repeated the computations 
for Rp = 400 (black curve) and Rp = 7000 (blue curve). The results presented in Fig. [3] 
can be summarized as follows. By fitting the three curves with their corresponding linear 
thermohaline diffusion coefficients (eq. [T5|) (dashed lines of the same colors), we determine 
very similar values of the effective finger aspect ratio, a ^ 0.4-0.5. This means that the 
effective salt-fingering /i-diffusivity can indeed be approximated by fllSp or, in other words, 
that Dp_/K oc l/Rp oc \V p\, at least within the investigated parameter range. Second, the 
value of a resulting from our direct 2D n umerical simulation s turn s out to be surprisingly 



close to the value of a = 0.5 advocated bv iKippenhahn et al.l (Il980l ). whereas it is an order 



of magnitude smaller than the value assumed by lUlrichl (119721 ). Third, the equilibrium value 
of Dp^ for our standard RGB case (i.e., for the values of relevant parameters in the vicinity 
of the HBS presented in Table [1]) is nearly an order of magnitude smaller than the value of 
Duax = 0.02i^ that satisfies the observational data (compare the red and green lines). Note 
that the latter was assumed to be constant during the upper RGB evolution, whereas Dp 
will obviously decline with time along with an evolutionary decrease of the envelope ^He 
abundance, hence |V^|, which will make it even more difficult for the ^He-driven thermo- 
haline convection to comply with observations (see the next section). Finally, the effective 
thermohaline heat diffusivity is found to be negligibly small compared to the radiative ther- 
mal diffusivity K in all the cases. Consequently, it will not influence the thermal structure 
of the radiative zone. 

The red curve from Fig. [3] is also plotted in the upper-left panel of Fig. HJ while the 
upper-right panel shows a snapshot of its corresponding salinity (/i in this case) fleld at 
equilibrium. A comparison of the latter panel with the upper panel of Fig. [2] leads us to the 
conclusion that, unlike the oceanic case, there are no visible vertically elongated structures 
(flngers) in the RGB case. This is the reason why we get a < 1 in the second case. The /x- 
fleld patterns in the upper-right panel of Fig. H] instead resemble disorganized turbulent-like 
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convection. A possible explanation of the difference in the salinity-field structures between 
the oceanic and RGB case can be given o n a basi s of res ults of the analysis of the equilibration 
of the salt-finger instability reported by iRadkd (|2010[ ). First of all, it should be noted that 
the difference between the two cases is evidently caused by their corresponding different 
values of the three dimensio nless pa r amete rs in equations (fT8l- [20|) : Pr, and r (the 



last three rows in Table [T]). iRadkd ( l2010l ) replaced the original non-linear equations by 



their weakly non-linear analogs using asymptotic expansions in which e = 1/Rp — t 1 
and then he solved them numerically for a wide range of Pr employing a computer code 
similar to that used by us. As a result, he has isolated two dominating mechanisms of salt- 
finger i ng eq uilibration. These are the triad interactions between various growing modes (e.g., 
Vallid I2OO6I ) and the adverse actio n of vertical s hear, spontaneously developing as a result 
of secondary salt-finger instability. iRadkd ( 120101 ) has come to the conclusion that, although 
both processes are essential, it is the shear instability that becomes critical at Pr <^ 1. 
Under these circumstances, the growth of fingering salinity and heat fluxes is limited by 
the energy dissipation sufficien t to balanc e the buoyancy forcing generated by the double- 
diffusive (primary) instability ( IShenI 119951 ). The fluxes generally decrease with the Prandtl 
number, although their dependence on Pr is non-monotonic. 

To verify these predictions for the RGB case, we have repeated our computations with 
the viscosity (Pr oc z/) artificially increased by factors of 10^ and 10^. The effective 11- 
diffusivities for these test cases, as well as their approximations by the linear diffusion coef- 
ficient f lT5]) . are plotted with blue and black solid and dashed curves in the upper-left panel 
of Fig. m We see that the increase of Pr indeed results in an increase of D^, the latter being 
proportional to the /i-fiux. The linear approximation of for ut = 10^ z/ = 4 x 10^ cm^ s~^ 
yields the effective finger aspect ratio a ~ 6.6 oq = 2.6, where = 0.4 for the standard 
RGB case. We use a subscript "t" in the expression for viscosity to emphasize the fact that 
this cannot be a microscopic viscosity anymore. Instead, it could be a turbulent viscosity 
like the one associated with rotation-driven instabilities. When u = 4 x 10^ cm^s~^, the 
/i-field patterns, like the value of a, appear to almost be identical to those obtained for the 
oceanic case (compare the lower-right panel of Fig. Hlwith upper panel of Fig. [2]). The case 
of z/t = 10^ u = A X 10^ cm^ s~^ corresponds to a transition from the jU-field structure at 
z/ = 4 X 10^ cm^ s~^ that looks more or less organized in the vertical direction to the struc- 
ture resembling a turbulent-like convection at z/ = 400 cm^ s~^. The most interesting pattern 
seen in the lower-left panel of Fig. H] is something like a vertically varying sinusoidal shear. 
This is most likely to be a manifestation of the mean fiel d effects (a secondary instability) 
that become critical at low Pr, according to iRadkd ( I2OIOI ). It is also interesting to compare 
our snapshot panels corresp ondin g to th e decreasing z/ with the sequence of T-field snap- 
shots presented in Fig. 2 by IShenI (jl995[ ). They look very similar. However, the important 
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difference is that Shen shows a time sequence along which the secondary instabihty, the one 
that leads to the equilibrium salt-fingering convection, is developing, whereas our sequence 
shows the equilibrium states already achieved for different values of Pr, at the highest of 
which the secondary instability appears to be strongly suppressed. 

Our results in the upper-left panel of Fig. HI in particular, the one represented by 
the black curve, could be used to support the hypothesis that the ^He-driven thermohaline 
convection is the main RGB extra-mixing me chanism, provided t hat a sufficiently strong 



turbulent viscosity, say of rotational origin (e.g.. lPalacios et al.ll2006l ). would be present in the 
radiative zones of RGB stars. However, such a speculation has to admit (and explain how it 
is possible) that the turbulence enhances only the viscosity but not the /x-diffusivity. Indeed, 
when we increase both v and v^^x by the same factor, the effective thermohaline /x-diffusivity 
returns close to its standard value (the green curve in the upper-left panel of Fig. H]). This is 
not surprising because an increase of z/moi (^s in the oceanic case) reduces the buoyancy force 
by making it faster for the difference in the chemical composition between rising and sinking 
blobs to be smoothed out horizontally (see the second term in the parentheses in eq. [T5|) . 
Furthermore, if we accept the hypothesis that turbulence in stellar radiative zones should be 
highly anisotropic, with its associat ed horizont al viscosity components strongly dominating 



over that in the vertical direction ( lZahnlll992l ). then the above speculation becomes even 



less likely, again provided that the chemical composition transport is accelerated by the 
turbulence proportionally to that of momentum. 



4. Post-Processing ID Simulations of the RGB Thermohaline Mixing 



The results of our 2D numerical simulations of thermohaline convection in the vicinity 
of the HBS in the bump luminosity RGB model predict the effective salt-fing er aspe c t rati o 
a = 0.4-0.5 that is an order of magnitu de smaller thaii the on e advocated by lUlrichI ( 119721 ). 
th e latter being needed t o repr oduce the iGratton et al.l (120001 ) observational data, according 
to ICharbonnel fc ZahnI ( 120071 ). The small value of a in the RGB case compared to its 



much higher value in the oceanic case is due to the more favourable conditions for the 
development of secondary instabilities, that strongly limit the growth of salt fingers, at low 
Prandtl numbers. 

However, our simulations are far from being perfect. Indeed, first of all, they are two- 
dimensional, in which case the oceanic salt-fingering fluxes were foun d to be underest i mate d 
by a factor of two to three compared to 3D numerical simulations ( iRadko fc Stern! Il999l ). 
For low Prandtl numbers corresponding to th e RGB case , differences between the 2D and 
3D solutions may be even more substantial (iRadkd 120101 ). Second, our simulations are 
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restricted to a small space domain surrounding the point of minimum Rp. Third, they 
do not take into account either nuclear reactions or a modification of Rp by thermohaline 
mixing, both of which may influence the growth of salt fingers. Finally, there still remains 
the possibility that the secondary instability associated with the mean field effects (vertical 
shear) can be suppressed in the RGB case, for example, by rot at ion- induced turbulence, 
such that the turbulent viscosity considerably exceeds the rate of turbulent mixing for some 
reason, which will result in a higher a. Therefore, we have decided to supplement our 2D 
numerical simulations of the RGB thermohaline convection with ID computations of the 
evolutionary changes of the surface chemical composition of RGB stars above the bump 
luminosity in which we model the ^He-driven thermohal ine convection using the li near- 



theory diffusion coeffi cient f|T5l) as in the work of, e.g. ICharbonnel &: ZahnI (120071 ) and 



Stancliffe et al.l (|2009[ ). However, given that the diffusion coefficient f[T^ is proportional to 



an extremely small quantity |V^| < 10~^, which is affected by the mixing itself, we decided 
that it would be very sensible to re-mesh the computational grid of full stellar evolution 
computations, and to perform our ID simulations on a fixed mesh in a post-processing way. 
Our goal is to see if we can adjust a value of a with which t he diffusion coe f ficien t (fT5|) 
will really be able to reproduce the ^^C/^^C, C, and N data of iGratton et all (120001 ). We 
want to do this also becau se th ere were some di s crepa ncies between the results reported by 



Charbonnel &: Zahij (120071) and IStancliffe et al.l (120091). on the one side, and those obtained 



and anticipated by iDenissenkov &: Pinsonneaultl (l2008bl ). on the other. 



4.1. Basic Equations and a Method of Their Solution 



Denissenkov fc Pinsonneaultl (j2008af ) have noticed that, in the absence of extra-mixing, 
the /i-profile in the radiative zone of an RGB star above the bump luminosity does not 
change very much if the mean molecular weight is plotted as a function of radius. The 
electron-degenerate He core in the center of the RGB star can be considered as a low-mass 
white dwarf, whose radius is known to weakly depend on its mass (R oc M^/^). The He core 
mass increases with time as the star climbs the RGB, thanks to the transformation of H 
into He taking place in the HBS atop the He core. This needs some fresh H to constantly 
be conveyed to the HBS from the base of convective envelope through the radiative zone. 
Note that, whereas the mass of the radiative zone is very small (~ 0.01-0.001 Mq), its 
radius extends from a value of the order of the Earth's radius at the He core boundary to 
a value of the order of the Sun's radius at the base of the convective envelope. A typical 
radial velocity of the mass infiow between the envelope and HBS, that feeds H to the HBS, 
is f ~ — (10~^-10~^) cms~^. Given these facts, it is natural to write and solve the nuclear 



-17- 



kinetics equations for the radiative zone in the Eulerian coordinates 

dVi , .dyi f dyi\ d f dyi 



dt'"'' dr~\dt J,,,i^ 5r dr ) ' ^^^^ 

where yi is the mole-per-gram abundance of the i"^ nuchde, is the thermohahne diffusion 
coefficient given by expression ( fT5|) with C = 2CKunzc = 27r^, and the first term on the 
right-hand side describes local changes of yi produced by nuclear reactions. 

Distributions of p(r), T(r), and other stellar structure parameters in the radiative zone 
necessary for calculations of r, D^, and the nuclear term in (1211) are taken from a couple of 
our MESA models separated by a sufficiently large distance in log L/Lq on the RGB, the first 
model being located immediately after the bump luminosity, where the RGB extra-mixing 
is supposed to commence. The models are used to interpolate the structure parameters in 
both r and t. The mass infiow rate can be estimated from the mass and energy conservation 
relations 

dr 1 dMr dMr L , , 

(22) 



dt Aur'^p dt ' dt ecN-^e' 



where is the envelope hydrogen mass fraction, ^cn ~ 6 x 10^*^ ergg~^ is the energy 
released per one gram of hydrogen burnt in the CN-branch of nuclear reactions, and L is 
the star's luminosity that is assumed to be constant in the radiative zone above the HBS. 
A comparison of these relations with real f-profiles from our RGB models has led to the 
following corrected relation 

that will be used in our ID simulations. 

Our test nucleosynthesis computations have shown that both the location and form of 
the /i depression are entirely determined by the reactions of ^He burning and those of the CN- 
branch of the CNO-cycle. Therefore, in our post-processing RGB extra-mixing simulations, 
we have included only the following reactions: ^He(^He,2p)^He, ^He(a, 7) ''Be(p,7)^B(/3"'"z/e) 
24He, i2c(p,7)i3N(/3+z/e)i3C, i3c(p,7)i4N, i^N(p,7)^^0(/ 3+^^J^^N, and ^/^N(p ,a)^^C. Their 



rates have been taken from the NACRE compilation by lAngulo et al.l (119991 ). The abun- 
dances of H, ^He, ^He, ^^C, ^^C, ^^N, and ^^N from the convective envelope of the MESA 
bump luminosity model are used as initial uniform conditions in the radiative zone. Simple 
boundary conditions, [dyi/dr) = 0, are applied at the bottom of the HBS, where X = 10^^, 
and at the star's surface. Mixing in the convective envelope is modeled with the diffusion 
coefficient -Dconv = 10^^ cm^ s~^, which keeps the envelope composition uniform. The system 
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of equations fl2Tl) has been solved numerically by a finite element method using the COMSOL 
Multiphysics software package. Initially, we substitute = u^^i into equations fl2T]) and 
solve them using the stellar structure parameters only from the first model for a long enough 
time to determine stationary abundance and fi profiles in the background of the radiative 
zone mass- inflow (e.g., panels a and c in Fig. |5]). After this, we switch on the ^He-driven 
thermohaline mixing and let it operate as our model evolves by adding the expression (fT5|) 
to Dfj^ and interpolating the relevant stellar structure parameters in time and radius between 
our two RGB models. Note that the addition of (ITS]) to the diffusion coefficient introduces 
an extra non-linearity in the PDEs fl2T|) because of the dependence 

^-°^ah^ = -^Z.(l + ^^)^^ah^' (24) 

i 

where Zi is the charge of the i^^ nuclide. It is this non-linearity associated with extremely 
small values of the term {din fi/dlnr) plus the frequent re-meshing of the radiative zone 
needed because the HBS constantly advances in mass that make the implementation of 
the ^He-driven thermohaline mixing difficult in full stellar evolution com put ations. The re 



searc hes who have done such calculations (e.g.. lCharbonnel &: Zahn]l2007l . and lStancliffe et al 



20091 ) did not provide details of their implementations making it impossible to try to repro- 
duce or assess their results. To eliminate a potential problem associated with the re-meshing, 
we solve the PDEs on a fixed mesh. To treat the non-linearity as precisely as possible, we 
use a large number of mesh points (approximately 1000) and set up very small tolerances in 
the time stepping algorithm of the COMSOL code. 



4.2. Solutions for the Low-Metallicity RGB Models 



A mean metallicity of upper RGB stars from the iGratton et al.l (120001 ) sample with 
[C/Fe] abundances clearly demonstrating an evolutionary decline corresponds to the heavy- 
element mass fraction Z ^ 0.0005, with a dispersion comparable to the mean value. Their 
masses as estimated using the Hipparcos parallaxes are centered around M ^ 0.85 Mq. 
Therefore, we will model the RGB extra-mixing in these stars using MESA RGB models 
that have Z = 0.001 and Z = 0.0001, and the masses M = 0.83 Mq and M = 0.85 M©, 
respectively. Panels a and c in Fig. |5]show the ^He, ^^C, ^^C, and n profiles in the vicinity 
of the HBS in the bump luminosity model with Z = 0.001 in the absence of extra-mixing. 
A rule of thumb for estimating a correct depth of the /i depression is the relation Afi ^ 
yU^AX(^He)/6, where AX('^He) can be replaced with the envelope ^He mass fraction with 
no harm done. Applying this rule for log^^o ^e(^He) ~ —2.8 and /i ~ fi^am = 0.595 (panels 
a and c), we obtain an estimate of Afi ^ —9 x 10^^ that is close to the maximum depth of 
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the n depression in panel c. 

The minimum of is reached at In rmin/i?© ~ —2.73, or rmmZ-R© ~ 0.065, where 
nin = ^(Atmin)- This is Very close t o the values reported in our previous publications (e.g.. 



Denissenkov fc Pinsonneaultl l2008bl ) . If we draw a vertical line through the point of /imin in 
panel c toward the panel a then we will come to our former conclusion that the ^He-driven 
thermohaline convection cannot explain the observed decline of [C/Fe] unless it penetrates 
below /imin- It turns out that such overshooting is indeed possible, provided that the thermo- 
haline mixing is much more efficient than the one predicted by our 2D numerical simulations. 
Panels b and d in Fig. [5] show the abundance and /i profiles that are obtained for a case of 
extra- mixing in which Z)^ has been modeled by equation f|T5|l with a = 5. Note that this 
finger aspect ratio is nearly ten times as large as the effective one estimated from Fig. [31 the 
resulting diffusion coefficient being almost two orders of magnitude higher. Panels b and d 
correspond to a short time after the mixing has been switched on. We see that the mixing 
extends down to the radius r^ix ~ .05 Rq. Again, this is close to the mixing depth used by 



Denissenkov fc Pinsonneaultl (j2008al ). Although the overshooting from rmin to rmix has been 
produced numerically, it has some physical justification, and can therefore be real. When 
the mixing is switched on, it will first steepen the /i-profile immediately to the right from 
'"min (panel c). Very soon, this will result in the formation of a discontinuity in the /^-profile 
at rmin, where material with a higher /i overlies material with a lower fi. Such a stratifica- 
tion is also subject to the double-diffusive instability, like the one with uniform salinity and 
temperature gradients considered previously. Hence, the material from the right will start 
to mix with the material from the left, thus pushing the /x-profile discontinuity to the left 
and, simultaneously, lowering it. This process will continue until the point rmix is reached 
at which the local reduction of fi by '^He burning, with ^He being conveyed by mixing from 
the envelope, is balanced by the increase of /i caused by the transformation of H into He in 
the CN-cycle. 

Fig. [6] compares the abundance and fi profiles for a = 5 (panels a and d), a = 6 (panels 
b and e), and a = 7 (panels c and f) in the same bump luminosity model. Black curves in 
the three upper panels show the ^^C profile in the model without mixing (the red curve in 
Fig. E^). One can see that the mixing with the higher finger aspect ratio penetrates deeper 
into the HBS and, at the same time, its diffusion coefficient is increased proportionally to 
(eq. [T5|) . However, this does not turn out to strongly change the evolutionary variations 
of the surface composition produced by the mixing with the large values of a. They are 
plotted with red (a = 5), green (a = 6), and blue (a = 7) curves in Fig. [TJ where crosses and 



triangles represent the iGratton et al.l ( l2000l ) data, the triangles showing upper limits to the 
i2q^i3q YdXio. For comparison, black curves in Fig. [7] correspond to the case of a = 1. Note 
that even this small finger aspect ratio still leads to a nearly four times larger thermohaline 
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diffusion coefficient than the one predicted by our 2D numerical simulations (the red curve 
in Fig. [3]). 



Denissenkov fc VandenBergj (120031 ) have demonstrated that an increase of the RGB 
extra-mixing rate by the factor of two should result in a strong enhancement of the evolu- 
tionary decline of the carbon abundance (see their Fig. 7). An increase of the mixing depth 
should also enhance the [C/Fe] depletion, though in a lesser proportion than that of the 
diffusion coefficient (their Fig. 8). In the present case, it turns out that, in spite of the fact 
that the increase of a from 5 to 7 doubles the diffusion coefficient (1151) and, at the same time, 
produces deeper mixing, this does not affect the [C/Fe] evolutionary decline (compare the 
red and blue curves in Fig. [7)d) very much. This difference is explained by a convergence 
of the diffusion coefficients f lT^ calculated for different but sufficiently large values of a to 
the same proffie (Fig. [H]), as ^He gets depleted in the envelope (Fig. [Tji). The faster and 
deeper mixing destroys ^He quicker than the slower and shallower mixing and, because the 
diffusion coefficient (fT5|) is indirectly proportional to the mass fraction of ^He left in the 
envelope (through the dependence A/i oc AX('^He)), the D ^-profile corresponding to the 
larger value of a quickly converges to that calculated for the smaller a. In other words, 
it can be said that the efficient ^He-driven thermohali ne mixing becomes self-quenching . 



On the contrary, the aforementioned results reported by iDenissenkov fc VandenBergl (j2003[ ) 



were obtained for constant r nixing rates. In connection t o this, it shou l d be noted that 



Denissenkov fc Pinsonneaultl (j2008bl ) have reproduced the iGratton et al.l (120001 ) data with 
the diffusion coefficient -Dmix = 0.02 K also assuming that it does not change with time. 
Therefore, given that the coefficient f lT^ rapidly decreases with time as a result of the self- 
quenching (Fig. [8]), its initial values have to be much larger than 0.02 (Fig. |8^). Hence, 
the green line in the upper panel of Fig. [3] with which we have compared our 2D numerical 
simulations of the thermohaline diffusion coefficient should be increased by approximately 
one order of magnitude. 

The floor of the depression is found to be almost flat in the case of Z = 0.0001 and 
M = 0.85 M0 (Fig. |9t). As a result, for the same value of a = 5, the mixing in this model 
penetrates deeper (with respect to the ^^C proffie in the unmixed model, shown with black 
curves in Figs. Ejo and [Hb) than in the model with the higher metallicity (compare panels b 
and d in the two figures). As in the previous case, the mixing depth increases with a (Fig. fTO]) . 
Given that the bump luminosity is higher for the lower metallicity model, the mixing in it 
starts with a larger initial K, hence with a higher oc K, because K oc L but, at the same 
time, the mixing has less time to accomplish its task. Evolutionary changes of the surface 
chemical composition produced by the mixing are shown in Fig. [TT] We see that the low 
[C/Fe] ratios in the most luminous RGB stars from the Gratton et al. sample are almost 
reproduced with the diffusion coefficient f|T5|) . provided that the finger aspect ratio can reach 
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a value of a > 7. Note, however, that this resuh is obtained for Z = 0.0001, whereas some 
of these stars have metalhcities closer to Z = 0.001, in which case the agreement between 
our model predictions and the observational [C/Fe] data is worse (Fig. Wp)- Nevertheless, 
given the model uncertainties, the approximations used in our ID simulations of the RGB 
thermohaline mixing, and the fact that the metalhcities of upper RGB stars from the Gratton 
et al. sample are distributed between Z f» 0.001 and a value of Z close to 0.0001, we find 
the general agreement between the linear theory with the high finger aspect ratios and 
observations to be satisfactory. Besides, the model confirms the observational inference that 
the effect of the RGB extra-mixing increases towards lower metalhcities. Thus, the results 
of our ID simulat ions of ^He-driven th e rmoh aline mixing in upper RGB stars go along with 
those reported by ICharbonnel &: ZahnI (120071 ). In particular, even our adjusted finger aspect 
ratio a = 7 results in a value of the total non-dimensional coefficient Ct = 27r^a^ = 967 that 
is very close to the value of Ct = 1000 used by them. 



4.3. Solutions for the Solar-Metallicity RGB Model 



The //-profile in the vicinity of the HBS in our MESA bump luminosity model with 
Z = 0.0188 and M = 1.5 Mq is plotted in Fig. [121 It is shown in the absence of mixing. It 
has two important differences, as compared to the //-profiles in the low-metallicity models, 
by being shallower and narrower. The first property means that the ^He-driven thermohaline 
mixing in the solar-met allicity RGB star should be much slower than that in low-metallicity 
red giants for the same finger aspect ratio. A comparison of the coefficients for our three 
MESA models calculated using equation (flSjl with a = 7 confirms this conclusion (Fig. [T2|) . 
The second property means that the mixing in the solar-metallicity model cannot penetrate 
as deep into the HBS as it did in the low-Z models. As a result, it cannot dredge up carbon 
depleted material and, therefore, does not affect the surface [C/Fe] abundance at all. The 
only visible abundance changes that it can cause are Li depletion, and a modest reduction of 
the carbon isotopic ratio (Fig. [T^ . It is interesting that the linear theory (eq. [TB]) with a = 7 
predicts that '^He-driven thermohaline mixing in the solar-metallicity RGB star should lead 
to ^^C/^^C ~ 12-14, which is very close to the abu ndance ratios measured in upper RGB 



stars in the open cluster M67 (IGilroy fc Brownlll991l). which has nearly the sol ar metallicity 



and an MS turn-off mass of about 1.3 (e.g., IVandenBerg fc StetsonI 12004 ) . This result 
shows that the ^He-driven thermohaline mixing modeled by the diffusion coefficient (IT^ 
with the finger aspect ratio a ^7 can indeed reproduce properly the evolutionary abundance 
ch anges observed in upp er RGB stars of all metalhcities, a conclusion similar to that obtained 
by ICharbonnel fc ZahnI ( 120071 ). Having said that, we emphasize that our 2D numerical 



simulations predict a < 1; in which case, the model of the RGB thermohaline mixing fails 
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to interpret the observations (the black curves in Figs. [6] and [TO 



Li-rich giants may pose another problem for the explanation of the RGB extra-mixing 
by the ^He-driven thermohaline convection. These stars have masses and metallicities close 
to those of our solar- metallicity MESA niodel, a nd most of them are loc a ted ne ar the bump 
luminosity ( ICharbonnel fc BalachandranI 120001 ). iDenissenkov fc Herwigl (120041 ) have shown 
that the anomalously la rge abundances of Li in t hese stars can be explained by the "^Be- 
transport" mechanism (jCameron fc Fowled Il97ll ) only if their radiative zones experience 
extra-mixing with a significantly enhanced diffusion coefficient (up to lO^^cm^s"^). The 
blue curve in Fig. [12] corresponds to a = 7 but it only reaches values of the order of ~ 
10®cm^s~^. To obtain the values necessary for efficient Li production, we have to assume 
a ~ 200, which does not seem realistic. Hence, we would rather attribute the enhanced 
extra-mixing in Li-rich giants to some alternative mixing process that replaces thermohaline 
convection in these stars. However, we think this would be inconsistent because it is difficult 
to understand why this alternative mixing process cannot be the universal one that operates 
both in the Li-rich giants, which represent a few percent of all upper RGB stars, and in all 
other u pper RGB stars in w hich its efficiency is reduced to D-^^i^ ~ 10^-10^ cm^ s~^ for some 



reason. 



Drake et al.l (120021 ) have reported that, in their sample of single K giants "among 
rapid {ysmi > 8kms~^) rotators, a very large proportion (~50%) are Li-rich giants" and 
that "this proportion is in contrast with a very low proportion (~2%) of Li-rich stars among 
the much more common slowly rotating K-giants". This correlation of the RGB extra- 
mi xing enhancement (needed for the Li enrichment) with the rapid rotation has been used 
by iDenissenkov fc Herwid (120041 ) to speculate that the RGB extra-mixing is actually driven 
by rotation and that the Li-rich gia nts had been spun u p as a result of their engulfing of 
massive planets. Note that, although jPalacios et al.l (120061 ) have claimed that the RGB extra- 
mixing cannot be associated with a pure rotational mechanism, there is still a po ssibility that 



an interaction of rotation and large-scale magnetic fie lds can drive the mixing (IBusso et al. 



200?! : iDenissenkov. Pinsonneault. fc MacGregorll2009l ). The only scenario that could explain 
this correlation in the model with thermohaline mixing would be to assume that the Li-rich 
giants swallowed low-mass MS companions that are enriched in ^He. This could explain the 
rapid rotation (a deposit of orbital angular momentum) and probably the speed-up of extra- 
mixing by the increased amount of ^He that was supplied externally. However, from our 
point of view, this scenario appears to be too complicated, and therefore highly improbable. 

Finally, a comparison of panels d in Figs. [6], [TOl and [12] leads to the conclusion that, 
whereas ^He gets strongly depleted in the low-metallicity RGB stars, its envelope abundance 
is reduced by a factor of only a few in the solar-metallicity model. On the one hand, this 
can be used as an argument against ^He-driven thermohaline convection as the mechanism 
of RGB extra-mixing if the presence of similar extra-mixing process in low-metallicity AGB 
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stars with masses M < 1 Mq surmised by iMasseron et al.l (120061 ) and iLebzelter et al.l (|2008[ ) 
will by confirmed by other observations. On the o ther hand, if the observational abun dance 
anomalies, in particular those in meteorites (e.g.. iNoUett. Busso. fc Wasserburgll2003l ). sug- 



gest the operation of an RGB-like ex tra-mixing only in the population I AGB stars with 
masses M > 1 Mq (IBusso et al.ll2010l ). then a more detailed analysis of the consequences of 
modest ^He depletion in the solar-metallicity RGB model for the following AGB thermoha- 
line mixing driven by ^H e burning is still needed, especiall y given the arguments against this 
hypothesis presented by iKarakas. Campbell, fc Stancliffd (120101 1 . We will provide such the 
analysis in our forthcoming paper. 



5. Discussion and Conclusions 



In this work, we have come to the following two conclusions that happen to contradict 
one another. On the one hand, the linear stability analysis of the Boussinesq equations 
([T]-[3]) with parameters set up to describe the growth of salt fingers driven by ^He burning 
in the vicinity of the HBS in a low-mass RGB star above the bump luminosity leads to the 
thermohaline diffusion coefficient (|T5|) . It models surprisingly well the observed evolutionary 
changes of the surface chemical composition in upper RGB stars of different metallicities, 
provided that we employ the empirically constrained finger as pect ratio Onhs ^. 7. In othe r 
words, we arrive at a solution very similar to that proposed by ICharbonnel fc ZahnI ( 120071 ). 
namely that the effects of the RGB extra-mixing can be reproduced in stellar evolutionary 
computations if the simple diffusion coefficient 



V,.ri - V, 



(25) 



ad 



where Oobs ^ 7, is used. This large value of a is close to the one advocated by lUlrichI (119721 ). 
On the other hand, our 2D numerical simulations of thermohaline convection for the same 
RGB parameter set (Table [1]) have shown that the effective finger aspect ratio in expression 
( |T5|) does not exceed a value of Opff ~ 0.5 . Inter estingly, this small value of a coincides with 
the one estimated by iKippenhahn et al.l ( Il980l ). Because of the dependence of on the 
square of a, the difference between the two diffusion coefficients calculated with a = Oobs 
and a = a^s exceeds two orders of magnitude. It is highly unlikely that our 2D numerical 
simulations have underestimated by this much. Therefore, we are inclined to conclude 
that RGB extra-mixing has nothing to do with salt-fingering transport and that alternative 
mixing mechanisms are worth investigating. 



There is a clear physical explanation as to why salt-fingering leads to the turbulent 
convection in the RGB case, while it produces vertically elongated quasi-laminar structures 
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Shen 


1995; 


Radko 


2010) 



different outcome is the large difference in the Prandtl number (Tabled]). The low Prandtl 
number (lower viscosity compared to heat diffusivity) in the RGB case favours the develop- 
ment of secondary instabilities, predominantly the one associated with the mean field effects 
of vertical shear, which do not allow the salt fingers to grow longer than their diameters in 
the vertical direction. We have demonstrated that the artificial increase of viscosity does 
stabilize the growth of salt fingers, as expected. In real RGB stars, such an increase could 
be associated with turbulence generated by rotation-driven instabilities, in which case the 
turbulent viscosity z/t may exceed the niicroscopic (molecular plus radiative) viscosity v by 



several orders of magnitude (e.g., IZahnlll992l : iPalacios et al.ll2006l ). However, it is difficult 



to imagine how the turbulence can affect only the viscosity without enhancing chemical 
mixing at the same time. Hence, we also have to replace v^^i by t't in the expression in 
the parenthesis in the thermohaline diffusion coefficient f[T^ that has been omitted in fl2S]) 
for the sake of simplicity, because it can be neglected in the absence of turbulence. This 
replacement counterbalances the stabilizing effect of higher turbulent viscosity b ecause tur- 
bulent mixing, especially when it prevails in the horizontal direction ( 1Zahnlll992l ). facilitates 
the smoothing out of the mean molecular weight contrast between the rising and sinking 
fingers, thus weakening the buoyancy force. In connection with this, it is important to note 
that the hypothesis that RGB extra-mixing is only partially executed by the ^He-driven 
thermohaline convection an d that some other naixing mechanisms of rotational or magnetic 
origin are assisting it (e.g., ICantiello fc Langerl 120101 ) . does not seem to be plausible. In- 
deed, we know from the comparison of simple mixing models with the observations that 
extra-mixing in the majori ty of low-mass upper RGB stars ne eds a diffusion coefficient of 
the order of Z^mix ~ 0.02 K (jPenissenkov fc Pinsonneaultll2008bl ). If this mixing is associated 
with a mechanism different from that of the salt-fingering convection then we have to use 
the ratio D^^^jK ^ 0.02 instead of v^oxjK in the parenthesis in the expression ( [T5|) . Given 
that Rp = (Vrad — Vad)/V^ > 1000 in the radiative zones of RGB stars, the expression in 
the parenthesis becomes negative, which means that the density stratification is now stable 
against the double-diffusive instability. Consequently, there are only two possibilities: RGB 
extra-mixing is either entirely executed by the ^He-driven thermohaline convection or it is 
the result of an entirely different mechanism. To answer this question with certainty, our 
next step is to carry out 3D numerical simulations of thermohaline convection for the RGB 
case. 
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Fig. 1. — Solid curves show growth rates for salt fingers of different diameters for the oceanic 
(upper panel) and RGB (lower panel) cases. Dotted, dot-dashed and dashed curves are 
reciprocals of microscopic viscous, thermal and mixing timescales for perturbations with the 
wave number k = 2iT/d. 
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Fig. 2. — A snapshot of our 1024 x 1024 numerical simulations of the oceanic salt-fingering 
at the equilibrium state (upper panel). A change of color from blue to red corresponds to an 
increase of salinity. Blue and red curves in the lower panel show the diffusion coefficients for 
the salinity and heat salt-fingering transport in units of the microscopic thermal diffusivity. 
The blue dashed line approximates the equihbrium value of Ds/kx- 
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Fig. 3. — Upper panel shows thermohaline diffusion coefficients (oscillating curves) and their 
corresponding equilibrium values (dashed lines) obtained in our 2D numerical simulations 
for the density ratios Rp = 400, 1700, and 7000, which are close to those found in the 
radiative zones of RGB stars (lower panel). The solid green line in the upper panel gives 
the empirically constrained rate of the RGB extra-mixing. The diffusion coefficients are 
divided by the thermal diffusivity. In the upper panel, the inset also shows the effective 
finger aspect ratios (a) that produce the corresponding equilibrium values (dashed lines) 
when C = 2CKunzc = Svr^ is substituted in the expression (|T5l) . 
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Fig. 4. — Snapshots and results of our 2D numerical simulations of thermohaline convection 
for the RGB case with the same density ratio Rp = 1700 but for different values of viscosity 
and mixing rate. Upper-right panel: the standard RGB case (see Table [1]). Lower- left and 
lower-right panels: the viscosity has been increased by the factors of 10^ and 10^ (assuming 
it is now of turbulent origin). Upper-left panel: the diffusion coefficients and effective finger 
aspect ratios corresponding to the cases presented in the other three panels. The green 
curve has the same parameter set as the blue curve except that its corresponding mixing 
rate (molecular diffusivity z/moi) has also been increased by the factor of 10^. 
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Fig. 5. — The '^He, ^^C, and ^^C abundance and fi profiles in the vicinity of the HBS in 
the bump luminosity model with Z = 0.001 and M = 0.83 M0 computed in the absence 
of extra- mixing (panels a and c). Panels b and d: the same profiles after they have been 
modified by RGB extra-mixing modeled by the diffusion coefficient (ITSl) with a = 5. The 
black curve in panel b is the same as the red curve from panel a. 
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Fig. 6. — The first pair of tlie upper and lower panels (a and d) is the same as panels b and 
d in the previous figure, while the other two pairs correspond to finger aspect ratios a = 6 
and a = 7, respectively. 
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Fig. 7. — The evolutionary changes of the surface ^^C/^'^C ratio (panel a) and the abundances 
of C (panel b) an d N (p anel c) in the field low-metallicity and low-mass RGB stars studied 
by iGratton et al.l (120001 ) (crosses and triangles, the latter showing upper limits). Curves 
represent theoretical reproductions by our model of RGB extra-mixing with Z = 0.001 and 
M = 0.83 Mq in which mixing is described by the diffusion coefficient (1151) with a = 1 (black 
curves), a = 2, 3, and 4 (yellow curves), a = 5 (red curves), a = 6 (green curves), and a = 7 
(blue curves). Panel d shows the unobserved decline of the envelope ^He abundances that is 
predicted theoretically. 
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Fig. 8. — Convergence of the diffusion coefficient profiles for the models presented by the 
same colors in the previous figure. Panels a, b, c, and d show how the results vary with time. 
Each panel is separated, in turn, by a time interval that is approximately equal to 20% of 
the total time for which the models have been evolved. 
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Fig. 9. — Same as in Fig. [5] but for the bump luminosity model with Z = 0.0001 and 
M = 0.85 Mo. 
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Fig. 10. — Same as in Fig. [6] but for the bump luminosity model with Z = 0.0001 and 
M = 0.85 Mo. 
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Fig. 11. — Same as in Fig. [TJ except that the yellow curves are omitted, but for the RGB 
model with Z = 0.0001 and M = 0.85 M,^. 
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Fig. 12. — The ii and profiles and the evolutionary changes of the ^^C/^'^C ratio and 
abundance of ^He computed for the RGB model with Z = 0.0188 and M = 1.5 using 
the diffusion coefficient ( lT5l) with a = 7 (red curves). The blue and green curves in panel 
b show the profiles from the bump luminosity models with Z = 0.0001 and Z = 0.001, 
respectively. 
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Table 1. Correspondence Between Salt-Fingering Parameters 







Oceanic Case 






RGB Case 




Parameter 


Notation Value (cgs) 




Notation 


Value (cgs) 


Viscosity 




10-2 


V 






4 X 10^ 


Thermal Diffusivity 


kr 


1.4 X 10-3 


K 






108 


Haline Diffusivity 


ks 


1.1 X 10-^ 








2 X 10^ 


Gravitational Acceleration 


9 


9.8 X 102 


9 






10^ 


Thermal Expansion 


a 


2 X lO^"* 


a « 






10-7 


Haline Contraction 




7.5 X 10-4 7oo"' 








1.7 


Temperature Gradient 


dTo 
dz 


3 X 10-3 


dTo 
dz 


(dT-- 
\ dz / 


) ad 


2 X 10-3 


Density Ratio 


Rp 


1.6 


Rp - 


= (v- 


Vad)/V^ 


1.7 X 103 


Prandtl Number {v/kx) 


Pr 


7 


Pr 






4 X 10-6 


Inverse Lewis Number (ks/kx) 


T 


8 X 10-3 


T 






2 X 10-6 



